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Abstract 

With the recent investigations of g*g* — r/W transition form factor and rj — rf mixing 
scheme, we present an updated study of the radiative decays J/^f — > in perturbative 
QCD. The decays are taken as a test ground for the g*g* — rf'' transition form factors 
and the 77 — rf mixing scheme. The form factors are found to be working for glunic rf^ 
productions and the mixing angle is constrained to be cf> = 35.1° ± 0.8°. 
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As it is well known that the heavy quarkonium decays to light mesons have played very 
important role in testing and understanding QCD from the very beginning. The decays J/^> — > 
77^7 are of great interests since they are closely related to the issues of rj' — rj mixing and 
g*g*—T)W transition form factors, which are very important ingredients for understanding many 
interesting hadronic phenomena of 77 and 77' productions. For example, it would be very useful 
for explaining the large branching ratio of strong penguin dominated decay B — > Krj'^ |2l Ell- 
in the literature, studies of the decays J/^f — > 77^7 are different each other from the 
treatments of formatting gluons to rf'\ namely, direct nonperturbative gg—f]^ coupling through 
strong anomalyjlj or two off-shell gluons coupled to 77W through quark loopj^]. In this letter, 
we will take the second approach which had been pioneered systematically within perturbative 
QCD by Korner, Kiihn, Krammer and Schneider(KKKS) |S] years ago. In Ref. [Sj, the non- 
relativistic quark model and the weak-binding approximation were used for both heavy and 
light mesons, and systematic helicity projectors were constructed to reduce loop integrations. 
In this work, we follow their approach. However, two improvements are included: 

• g*g* ~ couplings are improved to be relativistic transition form factors as advocated 
in Ref. [HI EJ in stead of non-relativistic modelling. 

• The 77' — 77 mixing scheme is also updated to the Feldmann-Kroll-Stech(FKS) mixing 
scheme |S|. 

In perturbative QCD approach, the decays are depicted by the Feynman diagrams in Fig.l. 
To calculate the amplitudes for the decays, we need to know how to deal with the dynamics 
of bound states. Generally, factorization are employed. Soft nonperturbative QCD bound 
state dynamics are factorized to the decay constants and the wave functions of J/\P and 77W 
which will convolute with the hard kernel induced the decay. We shall use the non-relativistic 
approximation for the heavy J/ 1 !/, but not for the light mesons 77' and 77. Although a rigorous 
theory from first principles for the light bound-states are still missing, some effective approaches 
are in progress. In recent years, it has been realized that a proper treatment of the 77 — 77' system 
requires a sharp distinction between the mixing states and the mixing properties of the decay 
constants^. Taking strange-nonstrange flavor basis for the 77 — 77' system and the mixing of the 
decay constants following the same pattern of the state mixing, FKS have found a dramatic 
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simplification. They also have tested their mixing scheme against experiment and determined 
corrections to the first order values of the basic parameters from phenomenology. 
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Figure 1: Lowest order QCD diagrams for J/^f — > 7/^7 decays. 

In FKS mixing scheme the parton Fock state decomposition can be expressed as 

I 77) = coscj) I r) q ) — sincp \ r) s ), 

I rf) = sincj) \ n q ) + coscj) \ rj s ) (1) 

where cj> is the mixing angle, \n q ) ~ f q $(x, fj)\uu + dd) / \/2 and |r/ s ) ~ f s ®{x, fJ>)\ss). The decay 
constants / g , / s and the mixing angle cj> have been constrained from the available experimental 
data, f q = (1.07 ± 0.02)/,., f s = (1.34 ± 0.06)/^, cp = 39.3° ± 1.0° jHj. 

Already in Ref . ^0] , Baier and Grozin have derived evolution equations for the distribution 
functions fi) to the first order of a s , which eigenfunctions are found to be 

&(x,n)=6x(l-x)ll+ ]T B n (p)C%(2x- 1) ) . (2) 

\ n=2,4,... / 

In the limit \x — > 00, the coefficients B n evolve to zero and $(x, fi) turns out to be <Pas — 
6x(l — x). When evolution equations run down to low energy scale, its quark contents mixed 
with glunic states. However, the gluon content enters the rf 1 ' wave function from next-to- 
leading order. This observation encourages the calculations of the g*g* — 77 ^ transition form 
factors similarly to the well known 7* — ir transition form factor at leading order, which read 

0E1E] 

Mp, = (g* a g* b \v {,) ) = -^a s 8 ah ie^ aj 3Q"Q p 2 F g * g *_ n{l) {Ql,Ql), 

F g*g*-V^ Qt) = 7^rU» I d X - n 2 1 7^ 2 ~- + (X^X), (3) 

y y ' 2N C ' Jo xQi + xQ 2 — xxm^ (l) + le 

where x = 1—x, f v ' = \plf q sin cj>+f s cos cj> and f v = V2f q cos cj>— f s sin 0. To the accuracy of this 
paper, cp^,) (x, fi) is taken to be the leading twist distribution functions(DAs) (x) = 6x(l—x) . 
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Using g*g* — rf'' in EqJ3]and following the procedure developed in Ref.jSl^^; ^ * s straight- 
forward to evaluate the amplitudes for the decays as depicted by Feynman diagrams in Fig.l. 
We get 

nV - rj'l) = ~ (IJ e 2 Q at(M v )a e ^&(l - z 2 )\H{z)\\ (4) 

where z = m v i/My, eQ is the heavy quark electric charge and | is the color factor. The 
dimensionless scalar function H(z) containing loop integrals is given by 

M v 1 1 f l AS f 4 k v k 2 (p-kq 2 - q-kq-p) 

H{Z) ~ 2p^kTQ^Jo dU ^' {U) J d ^^klkliukl + ukl-uum^y (5) 

where D\ = —k\-{k + k 2 ), D 2 = —k 2 -(k + ki), q — k± — k 2 , and p = k\ + k 2 . 

Obviously in EqEJ the k\ ■ k 2 numerator would cancel the n' form factor if it is taken be 
~ l/ki-k 2 , and the hard scattering kernel would not convolute with the distribution functions 

Of 7]' . 

With the help of the algebraic identities 
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q2 Ml + ml 



m 2 AD 1 + D 2 ) + M v (k 2 1 + k 2 2 ) , (6) 



q-p = kj + k 2 2 , kvh = \{p 2 - k\ - k\) = - l -{p-k + D 1 + D 2 ), (7) 

the integrand in H(z) can be decomposed into a sum of four, three and two-points functions 
which is presented in appendix A. In the calculation of the loop integrals, we have used dimen- 
sional regularization scheme and the methods developed in Ref.|12j 

For numerical results for the decays, we use T tot (J/fy) = (87 ± 5)i^ey|T3], fj/q, = 400Mev 
and a s (M J/$ ) = 0.2557 JH|. We get 

B th (J/m -> rfi) = 3.9 x 1(T 3 , (B exp (J/^> -> 7/7) = (4.3 ± 0.3) x 1(T 3 , PDG^) (8) 
B th (J/m -> 777) = 3.5 x 10~ 4 , (B exp ( J/V -> 777) = (8.6 ± 0.8) x 10^ 4 , PDG[T$) . (9) 

While B th (J/^> — > r) >r y) agrees with experiment, B th (J/fy — > 777) turns out to be too small. 
From the mixing scheme, it is easy to see that B th (J/^> — > 7/7) is insensitive to the mixing 
angle <fi when <p is about 35°, but B th (J/^ — > 7/7) is very sensitive to <fi. Take <fi = 35.3° fitted 
from 7/ — > p'-f and p — > 777 [H], we find 

B th (J/ty -> 77V) = 3.75 x 10~ 3 , B th (J/V ^777) = 7.3 x 10" 4 , (10) 



which agree with experimental results quite well. However, if we take a s (/i) = ot s m c , the results 
turn out to overshoot their experimental data. 

The most theoretical uncertainty may arise from the energy scale choice in a s (/x). Because 
our calculation is performed at the lowest order in QCD and there is no UV divergence in the 
loop diagram induced the decay, we don't have strong argument to choose a scale, as in usual 
case, to minimize the higher order corrections by setting logarithm to zero. Naively, the scale 
could be chosen from m c to mj/^. To reduce the scale dependence, we relate B{J/^ — > 77W7) 
to B(J/« -> ggg) 

BW * - „« 7 ) = T Mt^t± B(JI* ggg). (11) 
With the help of the known results ^H] 



Yiy^ggg) 10(tt 2 - 9) a 3 s (M) f 1 , a s (M) 



1 + 



we can get 



3 „ / _ , f2M\ 



-19.4 + -A) 1-16 + In 



1 V ' \MvJ 



(12) 



afTtJ M , 9 e£ a s (M)a e f 2 (l) (l - z 2 )\H{z)\ 



20^ - 9) M£ j + ^ L 19 4 + | A ^ lg + ln (|M 



(13) 



We will use the following relation and experimental data^2] for our numerical results 



B(J/ty —> ggg) = B(J/^ — > hadrons) — B(J/^5> — > virtual^ — > hadrons) 
= (0.877 + 0.005) - (0.17 + 0.02) 

= 0.707 + 0.025. (14) 

Taking = 35.3° and a s (M) = a s (m c ), we obtain 

B th (J/ty -> 7/7) = 4.17 x 10~ 3 , B th (J/V 777) = 8.16 x 10" 4 , (15) 

which agree with the experiment data. 

In Fig. 2, we display the ratio 7Zj/y = B(J/^/ — > r] /r y)/B(J/^/ — > 777) as a function of 0, in 
which we expect that the relativistic and the higher order QCD corrections may be concealed to 
large extent, so that the ratio could be predicted much more reliable than the two decay rates 
respectively. Comparing our results with the experimental measurement Rj/q, = 5.0 ± 0.6[T3*] 
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Figure 2: The ratio TZj/q, is shown by a solid curve as a function of (f> (in degree). The experimental 
data are shown by horizontal lines. The thicker solid horizontal line is its center value, thin horizonal 
dash lines are its error bars. 

as displayed by horizontal lines in Fig.2, we find = 35.1° ± 0.8° which is different from the 
value = 39.0° ± 1.6° by 2cr also determined from J/^> — > ??(?/)7!H] by using QCD anomaly 
dominance mechanism formula. To make clear the origin of the discrepancy between the two 
different determination of mixing angle <p, we recapitulate the key formula from the well known 
work of Novikov et al. [3] 

J/ * (0\GG\ v ) \pj ' 
This formula is frequently employed to determine 77 — 77' mixing angles in the literature. Techno- 
logically, the strong anomaly dominance is equivalent to the dominance of the ground state and 
the neglect of continuum contribution to dispersion relations as shown in Ref . jU E] • So far, 
considering the experimental and the theoretical uncertainties, the difference between the pre- 
dictions for IZj/y by the two mechanism is still marginal. Although we have improved g*g* — rf'^ 
couplings in Ref. jH| from non-relativistic to relativistic, there is still large room for theoreti- 
cal improvements which is very worth for further studying. We also note the CLEO/CESR-c 
project is going, where about one billion \1/ events would be produced. The refined measure- 
ments of these decays to be performed at CLEO-c will deepen our understanding of the two 
rf 1 ' production mechanisms. 

In this letter, we have studied the radiative decays J/ty —* 7/(77)7 in perturbative QCD. The 
relativistic g*g* — rf 1 ' transition form factors have been tested to be working for 77' production. 
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The mixing angle in FKS scheme is constrained to be = 35.1° ± 0.8°. This study encourages 
further applications of the form factor for rf'> production in hard processes. It is also very 
helpful for understanding the abnormal large rj' yields in B meson decays, which have caught 
many theoretical attentions|17| ITH] recently. 
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Appendix A 

In the evolution of the amplitudes for J/ty — ■> •yq^, we encounter the loop integral in Eq|^l 
which can be expanded in terms four, three and two points functions 



H(z) 



1 1 

16 1 - 2 s 



/ La+^V 4 ( m4jD oK z) - ^(1 - z 2 )m 4 v D b (u, *)) 

1 



v l - uz 2 ) A 3 C b (u, z)--(l- 2z 2 + uz 2 ) 4 3 C a (w, z) 
+ -U 2 (5 >, z) - B b (u, z) - B c (u, z) + B*(u, z) 



(17) 



with the following functions 

D a (u,z) = 



771y(l — U)UZ 2 (1 — Z 2 ) 

1 — uz 2 



Sp 1 



1 — uz 
u{l-z 2 ) i 



2m\n 1 



1 — uz 
u(l - z* 



— l€ 



+ Sp\l- 



2ni + In 



1 - (1 -2u)z 2 J 
( u-(\- 2u)z 2 ' 



1 — uz 2 



2ni In |1 - 
I I I ni + In 



1 — uz 



— ie 



Sp[l- 



u — (1 — 2u)z 2 



l-(l-2w)^ 2 J r \ u(l-z 2 ) 
f(l-z 2 )(u-(l-3u + 2u 2 )z 2 )\ ' 



uz z 



D b (u,z) 



4m^(l - (1 - 2u)z 2 )u{l - z 2 ) 



2Sp 



uz 1 



2Sp 



(l-2u)(l-z 2 Y 



u 



+ Sp 
+ In/ 



z 2 (u- (1 -3u + 2u 2 )z 2 ) J 

l-uz 2 A, ((l-z 2 ){u 



Sp - 



z 2 (u - (1 - 2u)z 2 



In 



(l-2u) 2 z 2 (l-z 2 ) \ 
u- (1- 3u + 2u 2 )z 2 ) 
l-3u + 2u 2 )z 2 ) . N 

n V 



UZ' 



(19) 



C b (u, z) 



CS(u,z) 



BS(u,z) 



B b (u,z) 
B c (u, z) 

B d (u,z) 



ITly J 



dy- 



ln 



Au 2 z 2 - y 2 (l - 2z 2 ) - 2uy(l - 3z 2 ) - ie 
( |/(2(1 - 2z 2 ) - y(l - 2z 2 ) - 2u(l - 3z 2 )) - ie\ 



my J 



dy- 



2(y{l-z 2 ) - 2v?z 2 ) -ie 
1 



In 

2 



y 2 + 2yu(l — z 2 ) + Auz 2 — ie 
y(y + 2u{\ -2z 2 ) + ie\ 



e 

-U 



Au 2 z 2 + ie J ' 

7s + In An + In /i 2 + 2 - In m v 

In (l - (1 - 2u)(l - 2UZ 2 



1 



(1 - 2u)(l -2wz 2 ) + ze 



- 7s + In 47r + In /i 2 + 2 - In m^, 



e 



7s + In 47r + In jj 2 + 2 - In m 2 / 



/.e 



(1 - 2u)(2uz 2 - 1) +ie 



In 



7s + In 4vr + In jj 2 + 2 - In mt - 

e 1 — 22^ 



1 - (1 - 2u){2uz 2 - 1) +ze 
2(1 -z 2 ^ 



(20) 



(21) 

(22) 
(23) 

(24) 



ln[2(l - z z )\ + ln(l - 2«(25) 



Where Sp(x) = Li 2 (x) is the Spence function and u — 1 
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